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Abstract 

The state evolution of the initial optical noon state is investigated. The residue entanglement of the state 
is calculated after it is damped by amplitude and phase damping. The relative entropy of entanglement of 
the damped state is exactly obtained. The performance of direct application of the damped noon state is 
compared with that of firstly distilling the docoherence damped state then applying it in measurement. 

1 Introduction 

Quantum entanglement between two or more particles has attracted great interest and produced many applica- 
tions in quantum information processing, such as quantum communication, quantum computation and quantum 
cryptography and quantum metrology. It has been known for some time that entangled states can be used to 
perform supersensitive measurements, for example in optical interferometry or atomic spectroscopy pQ [HI- 
The idea has been demonstrated for entangled states of two photons0, three photons and four photons |SJ. 
In the best case, the interferometric sensitivity can reach the quantum mechanical 'Heisenberg-limit' in entangle- 
ment enhanced measurement which overwhelms the classical shot noise limit. If <fi is the phase to be estimated, 
and N is the number of independent trials in the estimation, the classical shot noise limit is A(j) = 1/VN. 
When entangled state is used, the limit can be reduced to at most to 'Heisenberg-limit' A(j> = X/N, One of such 
entangled states is the so-called noon state 

\N::0) ab =^=(\N,0) ab + \0,N) ab ) (1) 

which describes two modes a, b in a superposition of distinct Fock states \n a — N,nt — 0) and \n a = 0, n& = N) . 
The applications of this state include quantum metrology and quantum lithography 8 . In all the 
applications of noon state, the decoherence of the state and the performance of the damped state are less 
concerned. Huelga et al considered the ion system in presence of decoherence 0. It is inevitable that quantum 
state interacts with environment which will cause the decoherence of the state, thus in quantum supersensitive 
measurements decoherence should be included. We in this paper will investigate the entanglement of optical 
noon state in presence of decoherence and the performance degradation. 



2 Decoherence 

A quantum state will undergo decoherence after preparation. The decoherece comes from the interaction with 
environment. For continuous variable (CV) system, two most popular decoherences are amplitude damping 
and phase damping. The master equation describing these two decoherences for the density operator p is |l()j 
[TT][T2] (in the interaction picture) 4? = (£i+£2)p,with C\ represents the amplitude damping concerning with 
vacuum environment, 

dp = ^2 y (2«ipat - a\a,p - pa\a,i), (2) 



1 



and £2 represents the phase damping, 

£-2P = X! ~f [2a\aipa\ a i ~ i a \ a i) 2 P ~ P{4, a if\ > ( 3 ) 

i 

with a$ the annihilation operation of i — th mode, and 1^ and 7$ are damping coefficients of i — th mode for 
amplitude and phase damping respectively. The solution to the master equation can be conveniently obtained 
by first transforming the equation to the diffusion equation of the characteristic function of the state, then solve 
the the differential equation of the characteristic function. The time evolution solution of he density operator 
can be recovered from the characteristic function. 

The characteristic function is defined as % = tr[pD(p)], where T>(p) — exp(/ia^ — p*a) is the displacement 
operator, with p = [pi,p2, ■ ■ ■ , Ps] ,a> — [^1,0,2, ■ ■ ■ ,a s ] T and the total number of modes is s. The diffusion 
eqation of the characteristic function will be ^Sj 

!~5£ r rfWlg|i + W"M + i£T.S. W 

3 3 3 

where pj — \pj \ e l6j . The solution is simply be 

X (p,P*,t) = fdx X (pe^ +lx ,p*e-^- lx ,0)Y[(2n lj tr 1 ^ (5) 
3 3 

exp[-^-i(l- e -^| W | 2 ]. 

with pe~^ +tx standing for {p,\e~~^ Hxi , p2£~^ htX2 i --- , /j, s e~~i~ +/lXa ). The time dependent state can be 

recovered by p = J[U 3 ^-]x(p, P*,t)D(-p) [H]. 
The characteristic function of the noon state is 

X(\N::0)) = (N :: 0| D(ji) \N :: 0) (6) 

= Ie"^[L w (| M i| 2 ) + L N (\p 2 \ 2 ) + ^((-P*iP2) N + (-P1P* 2 ) N % 

where Ln is the Laguerre polynomial of order N : Ln(z) = ELo^^C)- The time evolution of the 
characteristic function will be 

x(p,P*,t) = ~ e -^[MlMi| 2 e- rit )+MlM2|V r *') (7) 
where T = ^(Ti +T2), 7 = 3(71 +72)- The solution to the master equation of density operator will be 

N / \ 

ffi H [(1 - e-^) N - m e- mr ^ \m0) <m0| (8) 



2 L ^ \m 

ra=0 



+ (1 - e -T2t^N-m e -mr 2 t | 0to ^ 



+e- Nrt - N2 ^[\N0) (0N\ + \0N) (N0\}}. 
The integral on p is carried out by the technique of integral within ordered operators. 



3 The entanglement of the damped state 

The damped state p is a mixed state. According to Peres-Horodecki criterion, the state is always entangled. 
The entanglement of the state can be carried out if measured by relative entropy of entanglement. The relative 
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entropy of p with respect to a separable state a is S(p\\a) — Tr(p log 2 p — plog 2 <7), the relative entropy of 
entanglement of p is the minimization of S(p\\cr) over all separable state a. Let the extremal separable state 
that minimizes the relative entropy be a*. Denote p = coo 1 00) (00 1 + X)m=i( c ™o \mQ) (m0\ + co m \0m) (0m|) 
+c(\N0) (OAT | + | OAT) (AT0|). We will obtain a* by first suppose a* having a special form then prove that a* is 
extremal (see Appendix) . Suppose 



N 



1 00) (00 1 + ( d ™o \m0) (m0| + d 0m |0m) (0m|) 



(9) 



rn — 1 



-d(|AT0) (0AT| + \0N) (N0\) + d NN \NN) (NN\ . 



The function that should be minimized is 

- Tr(p log 2 a*) = -c 00 log 2 d 00 - ^ (c m0 log 2 d m0 + c m0 log 2 d m0 ) 



N-l 



(10) 



-Tr 



m— 1 

C N0 C 
c C 0N 



log 2 



dj\io d 
d doN 



The constraints are Tra* = 1 and d 2 — doo^zviv > the later specifies that the extremal separable state should 
be at the edge of the separable state set (e.g.^El)- For the general situation of cjvo 7^ cqn, the minimization 
problem have not an analytical solution 16 . When the amplitude damping is symmetric, we have Ti = T 2 = T 



(thus c m0 = c 0m ; m = 1, 



, N), the solution to the minimization problem is 
d m o = dom = Co m< for m = 1,-'- , N- 1; 



^00 

dozv 



c 2 c 00 



(cat + c 00 ) 2 - c 2 ' 
= rfyvo = cato — d?^N ; 



(cato +c 00 ) 2 - c 2 ' 

c(catq + Cqo)cqo 
(catq + c 00 ) 2 - c 2 ' 



The relative entropy of entanglement of state p is: E r (p) = Trp(\og 2 p — log 2 cr*) 
+c_ log with c± = catq ± c, d± = dA^o ± d. It can be written as 



coo log B 



(11) 
(12) 

(13) 



E r (p) = 2(coo + CAro)[l-i?2( 



Cqq + catq + c 
2(c o + cno) 



(14) 



where ^(e) — — elog 2 £ — (1 — e)log 2 (l — e) is the binary entropy function, and c o = (1 — e ) , cato = 
^e~ Nrt , c — ^e~ Nrt ~ N 7 *. It should be mentioned that the corresponding solution to the problem of two qubits 
system is known ^7] [T5] . 

Other entanglement measures are entanglement of formation and distillable entanglement. From the defini- 
tion of the entanglement of formation, it is easily to obtained an upper bound for the entanglement of formation, 
which is 



E f(p) = ( c no + c 0N )H 2 { 



1 + y/l- c 2 /{c m + cqn) 2 



(15) 



We suspect if this is just the entanglement of formation itself. For the symmetric amplitude damping, when 
there is not phase damping and N > 5, E^(p) is very close to E r (p). 

The distillable entanglement is lower bounded by the coherent information (hashing inequality) . The coherent 
information of the state is 



A' 



N 



I c (p) = -cno log 2 c N0 - c om) log 2 ( c 0m ) 

rn — 1 rn — 1 

+ ^2 Co m log 2 (c 0m ) + c + log 2 C+ + c_ log 2 C_ 



(16) 
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When only phase damping is considered, that is I"i = T 2 — 0, we have co m = for all m < N, cnq = 
5, c = \e~ N 7 *. The coherent information will be I c {p) = 1 — H 2 (^ + \ e ~ N 7 ')- Meanwhile the relative entropy 
of entanglement will also be E r (p) = 1 — H 2 {\ + \e~ N2lt ). Because distillable entanglement Ed(p) is upper 
bounded by the relative entropy of entanglement, we have I c (p) < Ed(p) < E r (p). Now I c (p) = E r (p), thus for 
the situation only phase damping, the distillable entanglement is: 

E d {p) = l-H 2 ( 1 -+ l -e- N ^). (17) 

4 Performance of damped state in entanglement enhanced phase 
measurement 

In the entanglement enhanced phase measurement, the measurement operator is A — \N0) (0N\ + \0N) (iVQ|0. 
We suppose the state undergo the phase shift just before the measurement apparatus. The state is mod- 
ified by the phase shift to p v = c 00 |00) (00| + £^ =1 (c m0 |m0) (m0| + c 0m |0m) (0m|) +c(e iN v \NQ) (QN\ + 
e~ tNip I OAT) (iV0|). It should be noted that the relative entropy of entanglement is not changed by the phase 



shiftHEI. Thus (A) = Trp v A = 2ccosN(p = e - Nrt ~ N f* cos Nip, AA = J (A 2 ) - (A) , the phase deviation is 

A A \/h(e- NTlt + e- Nr ^) - e - 2NTt ~ 2N2 ^ cos 2 Nip 

a<p = — — — = =— . (is) 

\d(A)/dip\ Ne -Nrt-N*?t jsinAT^I 

The minimal Aip will be at ip = (2k + l)^. Thus the bset phase measurement precision will be 



(A v ) best = I e ^yi(e^ lt + e Nr 2 ty (19) 

We may adopt another strategy to enhance the phase measurement precision with the mixed entangled state 
at hand. In this strategy, the damped state p is distilled to noon state, then we use the new noon state for 
measurement. The successful probability of distilling a noon state from the damped state is characterized by the 
distillable entanglement of p. Thus we use noon state in the measurement at a probability of p — Ed(p), the phase 
deviation is the Heisenberg limit (A<p) q = j^; In a probability of 1 — Ed(p) we have no quantum entanglement 
to enhance the measurement, the phase deviation will be the classical shot noise limit (Aip) c = -^=. The 

mixture of the two kind of measurements will have the phase deviation: (Aip)d = ^Jp(Aip) 2 + (1 — p)(Atp) 2 , 

(this may derived from the fact that the mixture probability distribution function fdijp) — pfqi^p) + (1 ~P)fc{'p), 
with f q ((p) and f c (<p) are probability distribution function of entanglement enhanced phase measurement and 
classical measurement, the two kinds of measurement have the same mean). In the distillation strategy, the 

phase deviation will be 

(Acp) d - ^Ed/N 2 + (1 - E d )/N. (20) 

We can compare the performances of the direct application of damped state to phase measurement and the 
distillation then measurement strategy. Let us firstly consider the situation of phase damping alone. We have 



^2 



[1 + (N — 1)H 2 {- + l e -^-*)] e - 2JV (21) 



( A( ^)Lt ' V 2 2 

As indicated in figure 1, the performance of direct application is better when the phase ip under measurement 
is near (2k + 1)tt/(2N). While distillation strategy is better when the phase ip is far from these values. 

In the amplitude damping situation, distillable entanglement is upper and lower bounded by the relative 
entropy of entanglement and the coherent information respectively. It is followed that the phase deviation is 
also upper and lower bounded. We have (Aip)di < (A(p)a < (A<p)d u with (Aip)di — \jE r jN 2 + (1 — E r )/N 
and (Aip) du — \/Ic/N 2 + (1 — I c )/N. In figure 2, we calculate the upper and lower bounds of resolution for the 
distillation strategy in situation of symmetric amplitude damping alone. We can see that the best resolution of 
the direct application of damped state is better than the distillation strategy. 
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0.2 0.4 0.6 0.8 1 

Distillable Entanglement 

Figure 1: The performance of state with phase damping alone in entanglement enhanced measurement. The 
phase resolution is defined as the inverse of phase deviation. The line groups from top to bottom are for N = 
4, 3, 2 respectively. In each group the solid lines from top to bottom are 1/Atp for ip = (1, 3/4, 1/2, l/4)7r/(2iV) 
respectively, the dashed line is for l/(Aip)d. 




1 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 

Relative entropy of Entanglement 

Figure 2: The performance of state with amplitude damping alone in entanglement enhanced measurement. 
The phase resolution is defined as the inverse of phase deviation. The line groups from top to bottom are for 
N = 4,3,2 respectively. In each group the solid line is for 1/ (Aip)i, est , the dotdash line is for l/(A<^)^,the 
dashed line is for l/{Aip)d u . 
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5 Conclusions 



The master equation of quantum continuous variable system is sovled in the case of simultaneous amplitude 
damping of vacuum environment and phase damping. When the initial state is a noon state, the exact expression 
of time dependent solution of density operator is obtained via the characteristic function method. An analytical 
formula is given for the relative entropy of entanglement of the damped state when the two modes of the 
noon state undergo the same amount of amplitude damping. In the asymmetric amplitude damping, the the 
relative entropy of entanglement can be calculated by numerically solving a group of algebraic equations. In 
the situation of phase damping alone, the exact distillable entanglement is given, which enables the comparison 
of two strategies of applying the damped state in phase measurement possible. When amplitude damping is 
present, we calculate the relative entropy of entanglement and coherent information of the damped state. We 
use these to specify the upper and lower bounds of the distillable entanglement. 

The performance of direct application of the damped state in phase measurement is better than that of 
firstly distilling the damped state then applying it in measurement when the phase under estimation is near 
(2k + l)n/(2N). While distillation strategy is better when the phase <p is far from these values. 
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Appendix: Proof of the extremal state 

We prove that a* is extremal by the fact that local minimum is also the global minimum when it is in regard 
to the relative entropy of entanglement^^. Hence we only need to prove that a* is the local minimal state. 
Let f(x,a*,<r) — S(p || (1 — x)a* + xa ) be the relative entropy of a state obtained by moving from a* towards 
some a. The derivative of / will be |19 | |18 | 

/>oo 

+ ty 1 p{a* + ty 1 8a)dt = TrBSa, 



g(o,.V) = £V 



where we denote (1 — x)a* +xa = a* — 8a, and the operator B has the following matrix elements in the eigenbasis 
{|Xn}}ofa*: 

B mn = (Xm\ B \Xn) = ° gXn _ ° SXm ( Xm \ p \ Xn ) . 

Xn Xm 

And when X m = Xn, the corresponding coefficient should be replaced with the limit value of Xn 1 ■ 
We should prove that for any separable state Xi 

df 

(0, er*, er) = TrASo- = TrA(cr* - a) > 0. 

where TrBo* = l[Ij3' But any a € V (separable state set) can be written in the form of a = Yl,iPi \ct l f} 1 ) (a l /3 l \ 
and so (0, a*, a) = X^iPifs (Oj 17 *; la^/?*) (a^/^l), The problem is reduced to prove that for any normalized 
pure state \ocj3) (a(3\ , 

(a/3\B\a/3) < 1. (Al) 

We here prove the situation of symmetric amplitude damping system, the more general proof for corresponding 
two qubits system was already been found 16 , and the proof of asymmetric amplitude damping system is similar 
to that of two qubits system. For symmetric amplitude damping , the operator 

N—i 

Coo 



B = 100) (001 + V (ImO) (m0\ + |0m) (0m|) 
«oo _ n 



m— 1 



+ + C T )(\m (N0\ + \QN) (0N\) + C T )(\NQ) (0JV| + \0N) (N0\). 
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Denote \a) = £ m a m \m) , thus (a/3| i? |a/3) = \a /3 \ + E*S(l«mA>l + |aoA*| ) + arXKftvl + 

|aAf/?o| 2 ) + 5(3^ - 2z)(&n PoUoPn + a (3 N a* N ^). Let K\ = 1 - Em=i |a m | 2 ,#2 = 1- Em=i IA™| 2 , (with 
< #i,# 2 < 0), and a = JK[ cos Q x , a w = V^i sin^e"' 1 , /3 = y/K^ cos 9 2 e t712 , (3 N = y/K^ sin 9 2 e^ . Denote 
rj = r/2 — r/i — r]3, after maximization on 77, we have 

(a(3\B\a(3) = ^-K x cos 2 6»i cos 2 6 2 + (1 - K 1 )K 2 cos 2 <9 2 + (1 - K 2 )K l cos 2 6>i 

"00 

+ ^#i# 2 sin 2 (0 1 + 2 ) + ^K 1 K 2 sin 2 (6 1 - 6 2 ). 

Suppose the extremal value is achieved at some K\ 7^ 0, 1, by derivative on K\, we obtain the supposed (may 
not exist) extremal value (a(3\ B |a/3) = K 2 cos 2 6 2 < 1. What left is to verify that when Ki,K 2 = or 1, 
(a/3 1 B I a/3) will not exceeds 1. The nontrivial situation is K\ = K 2 = 1, thus 

(a/3| B I a/3) = -^-(cos^i + cos0 2 ) 2 + 7^7- sin 2 </> x + —-^ sin 2 </> 2 , 
where </>i : 2 = #1 ± 6> 2 . By using the fact that ^(^7 + 3—) = 2 ^+^~ , We can obtain the maximum value as 

(a/3| B |a/3) ro = i(^ + ^)=l. 
Hence inequality (|A1(I is proved. So that tr* is the extremal state that minimizes the relative entropy. 
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